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Abstract 

We have calculated the 6s — 7s parity nonconserving (PNC) El transition 
amplitude, Ep^c, in cesium. We have used an improved all-order tech- 
nique in the calculation of the correlations and have included all signifi- 
cant contributions to Epj^c- Our final value Ep^c = 0.904^1 ± 0.5%^ x 
10 _11 ieaB(— Qw/N) has half the uncertainty claimed in old calculations used 
for the interpretation of Cs PNC experiments. The resulting nuclear weak 
charge Qw for Cs deviates by about 2a from the value predicted by the stan- 
dard model. 
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I. INTRODUCTION 



There is an ongoing discussion on whether measurements of parity nonconservation 
(PNC) in atoms can be used to study new physics beyond the standard model (for a history 
of PNC in atoms, see, e.g., the book [|TJ or review 0). Such a possibility relies on the accu- 
racy of the measurements and the theoretical analysis. The best accuracy so far has been 
achieved for the cesium atom. In 1997, the Boulder group measured the PNC amplitude 
in cesium to an accuracy of 0.35% p|. The best calculations available at that time were 
published by the Novosibirsk group in 1989 [|J and the Notre Dame group in 1990 ||. Both 
works claimed an accuracy of 1%. This claim of 1% accuracy was based in part on the com- 
parison of the experimental and theoretical values relevant to the PNC amplitude (energies, 
electromagnetic El transition amplitudes, hyperfme structure). Later, new, more accurate 
measurements (see, e.g. [§,0) led to better agreement between theory and experiment for 
the El transition amplitudes. This allowed Bennett and Wieman to suggest that the actual 
accuracy of the PNC calculations in cesium is 0.4% ||. With this accuracy, the value of the 
nuclear weak charge Qw of the cesium atom which follows from the measurements and 
calculations [|],[5] deviates from the standard model value by 2.5a ||. The implications of 
this deviation for physics beyond the standard model have been examined in several works 

This result generated many recent works revisiting calculations of PNC in cesium. Dere- 
vianko |T2|] demonstrated that the contribution of the Breit interaction to the PNC amplitude 
Epnc (-0.6%) is substantially larger than previous estimates and reduces the deviation from 
the standard model. The Breit contribution was neglected in [|J and underestimated in || . 
The result of Derevianko was later confirmed in independent calculations [13,14] . A new 



many-body calculation of Ep^c was performed by Kozlov et al. [Q. It was calculated in 
second-order in the residual interaction with averaged screening factors. The result is in 
excellent agreement with similar calculations by Blundell et al. || and differs by less than 
0.5% from the more complete ("all-orders") calculations ||{|]. 

Radiative corrections to the weak charge of order a were calculated a long time ago 
in [15] (see also [Pj[). However, there are important "strong field" radiative corrections of 



order Za 2 and Z 2 a 3 \ia (\/R n ) which have recently been considered in JIT] (note that the 
latter correction is larger). Here A is the electron Compton wavelength, and R is the nuclear 
radius. The main contribution to Ep^c was found to be about 0.4%. This contribution 
originates from the radiative corrections to the weak matrix element due to the Uehling 
potential (recently this contribution was calculated numerically in |T8||). 

Among other corrections considered one should mention a correction for the neutron 
distribution |T9l which is quite small: -0.2% of E PNC - 



Note that Breit, Uehling, and neutron skin corrections are smaller than the uncertainty 
of 1% claimed for the calculated Ep^c amplitude in works @J§[. This uncertainty was 
mostly associated with the accuracy of the calculations of correlation corrections. Therefore, 
there is a very important question as to whether calculations of correlations were, or can be, 
performed to better than 1% accuracy. The analysis of the accuracy involves a comparison of 
calculated and measured atomic quantities and some additional "internal" tests, such as, e.g., 
checking the stability of the results against variation of certain parameters. Comparison with 
experiment alone is an important but not sufficient part of the analysis. This is especially 
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true when better accuracy is claimed for calculations performed many years ago and all 
the details can hardly be reconstructed. In our view, the only reliable way to improve 
the accuracy is to do the calculations again, trying to improve the method and numerical 
procedures at every stage, and repeat the analysis of the accuracy That is what we do in 
the present paper. 

We have made several important improvements to the method developed in 1989 [|J. 
First, we calculated a new series of higher-order correlation diagrams which account for the 
effect of screening of the Coulomb interaction in the exchange correlation diagram. In our 
earlier work [|J the screening effect was calculated for the direct correlation diagram only, 
while screening factors were used for the exchange diagram. The values of these factors 
were found by looking at the effect of screening in the direct diagram; however, the effect of 
screening in the exchange diagram may differ. In our present work we calculate the effect of 
screening for both diagrams. 

Second, we have made significant improvements to the numerical procedures at every 
stage of the calculations. This involves, for instance, using a more dense coordinate grid, 
including more terms into the summation over virtual states, and restoring the lower Dirac 
component of the wave function everywhere. We have also used a more accurate nuclear 
charge distribution. 

Our result for the PNC amplitude in cesium (without Breit, radiative, and neutron 
distribution corrections) Ep^c = 0.908^1 ± 0.5%) x lCT^ea^— Qw/N) coincides exactly 
with the old result M but has a much smaller uncertainty. (N is the neutron number.) 
To avoid misunderstanding, we should note that particular contributions to this value are 
slightly different in the old and the new calculations (for example, the Hartree-Fock value of 
the amplitude has changed by 0.3%). Therefore, we are indeed talking about a new result. 

We have also calculated the contributions of the Breit interaction (-0.61%) and the 
Uehling potential (0.40%) to Ep^c- There is very good agreement between different 
calculations for these values and they can be considered well-established. The resulting 
value for the PNC amplitude is E PNC = 0.906 (l ± 0.5%) x W~ n iea B (-Q w /N). The 
neutron-distribution correction, estimated as (—0.2 ± 0.1)%, reduces our final value to 
E PNC = 0.904(l ± 0.5%) x 10- n iea B (-Q w /N). 

II. METHOD OF CALCULATION 

The method we have used here is very similar to that used in our 1989 work ||. As 
a zero approximation we use the relativistic Hartree-Fock method. The perturbation in 
our approach is the residual interaction - the difference between the exact and Hartree- 
Fock Hamiltonians. For summing of diagrams we use a combination of the "correlation 
potential method" , which is a way to treat correlations by introducing a single-electron 
operator (correlation potential) E, and "perturbation theory in the screened electron residual 
interaction" |24|,|25]], which is used to calculate E. The method is an all-order technique, in 
terms of treating correlations, and is not a version of the popular coupled-cluster approach. 
The dominating sequences of higher-order correlation diagrams correspond to real physical 
phenomena like collective screening of the Coulomb interaction between electrons and the 
hole-particle interaction. They are included in all orders in our technique. An important 
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strong point of the method is that its complexity does not go beyond the calculation of the 
correlation potential E. When E is ready, the calculation of energies and matrix elements is 
relatively simple. E is used to calculate single-electron Brueckner orbitals. The calculation 
of matrix elements with Brueckner orbitals already includes most of the correlations. There 
are additional "non-Brueckner" contributions to matrix elements such as structural radiation 
and normalization of many-electron states. These contributions are small for s and p states 
and in most cases they can be expressed in terms of derivatives of E. This makes the 
calculation of hyperfme structure, transition amplitudes, etc. similar to the calculation of 
energies. 

The method works very well for alkaline atoms. The calculated energy levels and transi- 
tion amplitudes deviate by a fraction of a percent from experiment. The hyperfme structure 
(hfs) of s and p states has an accuracy of better than 1%. 

A. Correlation potential method 

The Cs atom has one electron above closed shells. Therefore, it is natural to start the 
calculations from the relativistic Hartree-Fock (RHF) method in the V^ -1 approximation 
(one electron in the field of the N — 1 core electrons). The single-electron RHF Hamiltonian 
is 

H = ca-p+((3- l)c 2 - Za/r + V N ^ , (1) 

a and (5 are Dirac matrices and p is the electron momentum. For cesium the accuracy of 
the RHF energies is of the order of 10%. 

To improve the accuracy one needs to include correlations. In order to do so we introduce 
a "correlation potential" E [^3] . E is a single-electron energy-dependent non-local operator 
defined in such a way that its average value over a single-electron state a is the correlation 
correction to the energy of this state, 

Ae a ={a\E(e a )\a). (2) 

The calculation of E will be described in the next section. E is added to the Hartree-Fock 
Hamiltonian to calculate single-electron states of the external electron, 

(H Q + E - e)^ Br = 0. (3) 

There are two major benefits from inclusion of E in the equations. First, the iteration of E 
is an important higher-order effect which gives a significant contribution to the energy and 
wave function. Second, by solving Eq. ([|) we obtain single-electron orbitals which are called 
Brueckner orbitals which already include correlations. These orbitals can then be used to 
calculate matrix elements for hfs, PNC, etc. 

When appropriate higher-order diagrams are included into the calculation of E, as will 
be described in the next section, by solving Eq. (|3]) we obtain energies of s and p states of 
alkaline atoms with an accuracy of about 0.1%. This is a radical improvement over the 10% 
accuracy of the Hartree-Fock approximation. 

To accurately calculate the interaction of external fields with atomic electrons one needs 
also to take into account the core polarization effect. We do this using the time-dependent 
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Hartree-Fock (TDHF) method (which is equivalent to the random-phase approximation with 
exchange). In this method every single-electron orbital has the form 

where ip n is the unperturbed Hartree-Fock orbital (or Brueckner one for states above the 
core) and 6ip n is a correction to the orbital caused by the external field. To find 5ip n one 
needs to solve the equations 

(H - e)8^ n = -(8e + h + 8V h )^ n . (4) 

Here h is the operator of the external field, SVh is the correction to the Hartree-Fock potential 
due to modification of the core states in the external field, and 5e is the correction to the 
energy. 

Equations (f|) are first solved self-consistently for core states to obtain SVh- Then Brueck- 
ner orbitals are used to calculate matrix elements between states of the external electron 

M ab = (^\h + SV h \^}- (5) 

This expression includes most of the correlations and already gives hfs and transition am- 
plitudes in alkaline atoms to an accuracy close to 1%. There are also contributions arising 
from the normalization of states and structural radiation P3] . When these contributions are 
included the accuracy improves. 



B. Calculation of the correlation potential 

We use the Feynman diagram technique to calculate the correlation potential S. The 
main reason for this is that this technique is more convenient for inclusion of higher-order 
diagrams. The drawback of it is the necessity to integrate over frequencies numerically. 

In the lowest (second) order in the residual Coulomb interaction, £ is given by the direct 
and exchange diagrams in Fig. [|. The corresponding mathematical expressions are 

t^\e,ri,rj) = [ ^Gij(e + u)Q im U mn (uj)Q nj , (6) 
E^(e, r f , r ) = V / f ^-^4i»G jm (e + ^i)G mn (e + u x + u 2 )Q mj G nj (e + u 2 ), (7) 

where summation over m and n is a numerical implementation of the integration over r m and 
r n , Q is the Coulomb interaction, G is the Hartree-Fock Green function, IT is the polarization 
operator: 

Qij = (8) 

^(e) = E 7 ^ + LSa, ( 9 ) 
flijiu) =/gG y (w + a)Gtf(oO =En^ n (G i ,(e n + uj) + G l3 (e n -u J ))^ n . (10) 

Here |n) are core RHF states and (7) are RHF excited states; S shows the pole passing rule for 
integration over frequencies. (Let us remind the reader that we use the Feynman diagram 
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technique to treat a basically non-relativistic correlation problem, e.g., our polarization 
operator ( fL0|) represents the polarization of the atomic core, not the vacuum.) 

The integration over frequencies in expressions (^j) and (|7|) can easily be performed an- 
alytically, reducing the diagrams in Fig. [I] to the four usual Goldstone diagrams presented, 
e.g., in work ||23|| . When higher-orders are included in £, only the integration for the polar- 
ization operator (pi]) can still be done analytically. The other integrals have to be calculated 
numerically. 

As was demonstrated in our earlier work p|J , the most important higher-order correlation 



contributions to S come from screening of the Coulomb interaction between atomic electrons 
by other electrons and from the hole-particle interaction in the polarization operator (another 
important higher-order effect is the iteration of E, which is included separately, by inserting 
£ into equations for single-electron orbitals; see previous section). 

The hole-particle interaction in the polarization operator II (Fig. 0) is included by re- 
placing the V^ -1 Hartree-Fock potential in the calculation of II by the V N ~ 2 potential. 
Since the calculation involves iterations of the RHF equation, the hole-particle interaction 
is included in all orders. 

The screening of the Coulomb interaction is a more complicated effect. The correspond- 
ing higher-order diagrams can be obtained from the second-order diagrams (Fig. [TJ) by the 
insertion of a series of hole-particle loops into every Coulomb line. Therefore, instead of, 
e.g., one direct diagram (Fig. [I]) we have an infinite chain of higher-order diagrams (Fig. 
One can see that an internal part of this chain (before the integration over u) forms a matrix 
geometric progression 

fi + ngfi + fiQnQn + . . . 

The sum of this progression is 

= h(u)[l-Qfl(tu)]-\ (11) 

and can be called the "screened polarization operator" . Screening of the Coulomb interac- 
tion in the direct correlation diagram is included by replacing the unscreened polarization 
operator II in Eq. (||) by the screened polarization operator tt: 

£ d (e, r h Tj) =J2[ 7T- G v( e + u)QimKmn(v)Qnj- ( 12 ) 

It is convenient also to introduce an operator of the screened Coulomb interaction 

Q(u) =Q(l-Qn(co))- 1 . (13) 

This expression can be obtained in a similar way to the screened polarization operator 7r by 
summing a matrix geometric progression corresponding to the chain of diagrams in Fig. |j. 
An expression for the direct correlation diagram can be re- written in terms of Q: 

Sd(e, n, Tj) =^2 [ 7^-Gij(e + uj)Q irn tl mn (u)Q nj (uj). (14) 

nm J ^ 

The screened Coulomb interaction operator Q can also be used to include screening into the 
exchange correlation diagram. This is done by substituting Q — > Q(uS) in equation (0) 
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^ex( e ; r ii r j) — 

yy r r dui d ^ 2 Q. n ^Q. m ^ + uj^Gmnie + u)\ + u; 2 )Q m j(^2)G nj (e + u 2 ). (15) 
mn J J 2ti 2tt 

The second-order correlation potential with the hole-particle interaction and screening 
of the Coulomb interaction included in all orders is depicted in Fig. ||. 

In our 1989 work [|J we used an approximate expression for the screened Coulomb inter- 
action in the exchange correlation diagram 

Qk ~ fkQk, (16) 

where k is the multipolarity of the Coulomb interaction and fk is the screening factor. In 
this expression the dependence of the screening on frequency is neglected. The values of the 
screening factors were obtained by calculating the direct diagram with and without screening. 
The use of average screening factors which do not depend on frequency significantly simplified 
the calculation of the exchange diagram. Like in pure second-order, the Goldstone diagram 
technique and direct summation over a complete set of single-electron states were used. No 
integration over frequencies was needed. 

In the present work we do the full-scale calculation of the exchange diagram using 
Eq. (0). This makes the method totally ab initio since no screening factors are used. 
The drawback of the method is the need to do double integration over frequencies numer- 
ically (note that there is only a single numerical integration over frequencies in the direct 
diagram) . 



C. PNC amplitude 

To calculate the PNC 6s — 7s amplitude in Cs one needs to include two external fields 
acting on the atomic electrons: the weak field of the nucleus and the electric dipole (El) 
field of the photon. The nuclear spin-independent weak interaction of an electron with the 
nucleus is 

G 

Hw = T^p(r)Qwl5 (17) 

where Gp is the Fermi constant, Qw is the weak charge of the nucleus, 75 is a Dirac matrix, 
and p(r) is the nuclear density. The El Hamiltonian is 

H E1 = -d-E{e- lu)t + e iu}t ), (18) 

where d is the dipole moment operator. 

In the TDHF method, a single-electron wave function is 

if, = ^ + 5i/j + Xe~ lwt + Ye iuJt + 5Xe~ iuJt + 5Ye iuJt , (19) 

where 5ip is the correction due to the weak interaction acting alone, X and Y are corrections 
due to the photon field acting alone, and 5X and 5Y are corrections due to both fields acting 
simultaneously. These corrections are found by solving self-consistently the system of the 
TDHF equations for the core states 
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(H - e)5^ = -(H w + 5V W )^ } (20) 

(H -e- u)X = -{Hex + 6V E1 )^, (21) 

(H - e + u)Y = -(H E1 + 5V E1 )iP, (22) 

(H -e- uj)5X = -5V E1 5tfj - 5V W X - 5V EIW ^, (23) 

(H - e + u)SY = -5V El 5^ - 5V^Y - 5V E1W ^, (24) 

where SVw and 8V E \ are corrections to the core potential due to the weak and El interactions, 
respectively, and 8V E xw is the correction to the core potential due to the simultaneous action 
of the weak field and the electric field of the photon. 

The TDHF contribution to Ep^c between the states 6s and 7s is given by 

E T pnc F = (MHm + 5V E1 \5ip Gs ) + (tP 7s \H w + 6V w \X 6s ) + (ip 7s \8V Elw \ip 6s ) • (25) 

The corrections 8tpQ S and X 6s are found by solving the equations ( PU| - f2"l"] ) in the field of the 
frozen core (of course, the amplitude ( |25| ) can instead be expressed in terms of corrections 
to ip 7s ). 

If we use Brueckner orbitals instead of RHF orbitals to calculate the PNC amplitude 
in Eq. (|25| ) we include all-orders in £ contributions to Ep^c- However, the correlation 
potential is energy-dependent, S = S(e), which means that S operators for the 6s and 7s 
states are different. We should consider the proper energy-dependence at least in first-order 
in E (higher-order corrections are small and the proper energy- dependence is not important 
for them). The first-order in E corrections to Ep^c are presented diagrammatically in Fig. |6|. 
We can write these as 

{^ 7s \t s (e 7s )\5X 6s ) + (5^ 7s \t p (e 7s )\X 6s ) + (5Y 7s \t s (e 6s )\^ 6s ) + (Y 7s \t p (e 6s )\5^ 6s ) . (26) 

The non-linear in E contribution can be found by subtracting from the all-orders result the 
first-order value found in the same method. 

The correlation corrections to Ep^c we have considered so far are usually called 
"Brueckner-type" corrections. (In this case the external field interacts with the external 
electron lines.) There are also contributions to Ep^c i n which the external field acts inside 
the correlation potential (see Fig. |7|). Those diagrams in which the El interaction occurs 
in the internal lines are known as "structural radiation", while those in which the weak 
interaction occurs in the internal lines are known as the "weak correlation potential" . There 
is another second-order correction to the amplitudes which arises from the normalization of 
states [23|. The structural radiation, weak correlation potential, and normalization contri- 



butions are suppressed by the small parameter E ext /E coie ~ 1/10, where E ext and E int are 
excitation energies of the external and core electrons, respectively. 

We have also included a correction to Ep^c due to the Breit interaction. We calculated 
this in a way similar to the earlier work The main difference is that in the present 

calculations we have also included the Breit contribution to the last term in Eq. (E5|). This 
makes the calculations more consistent but doesn't change the result significantly. 

We postpone the analysis of radiative corrections until Section [TV]. 
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III. 65 - 7S PNC AMPLITUDE 



The results of our calculation for the 6s — 7s PNC amplitude are presented in Table |. 
Notice that the time-dependent Hartree-Fock value gives a contribution to the total am- 
plitude of about 98%. The point is that there is a strong cancellation of the correlation 
corrections to the PNC amplitude. The stability of the PNC amplitude compared to other 
quantities in which the correlation corrections are large will be discussed in more detail in 
Section |V[ Notice that the values do not differ significantly from our 1989 results (compare 
the 1989 final result 0.908 x lQ- u iea B (-Q w /N) with "Subtotal" of Table | for the current 
calculation). The new series of higher-order diagrams and higher numerical accuracy of the 
current work has therefore not changed the previous result (note, however, that particular 
contributions are slightly different). 

The mixed-states approach has also been performed in pfl and to determine the 



PNC amplitude in cesium. However, in these works the screening of the electron-electron 
interaction was included in a simplified way. In |J empirical screening factors were placed 
before the second-order correlation corrections E'- 2 *' to fit the experimental values of energies. 
Kozlov et al. WM introduced screening factors based on average screening factors calculated 



for the Coulomb integrals between valence electron states. The results obtained by these 
groups (without the Breit interaction, i.e., corresponding to the Subtotal of Table |) are 
0.904 H and 0.905 |14|. To be sure that we understand the difference between these values 



and our value, we performed a pure second-order (i.e., using T,^) calculation and fitted the 
energies (as was done in |J) and reproduced their result, 0.904. 

In the work || a calculation using the sum-over-states method was also performed. In 
the sum-over- states approach the 6s — 7s PNC amplitude is expressed in the form 



v ^ A7s\H w \np) (np\H E1 \Qs) (7s\H E1 \np}(np\H w \6s) \ 
pnc 2^V Ei - E Ea - E > ' 

^"s J-'np J-'bs J -'np 



The authors of reference H include single, double, and selected triple excitations into their 
wave functions. Note, however, that even if wave functions of 6s, 7s, and intermediate np 
states are calculated exactly (i.e., with all configuration mixing included) there are still some 
missed contributions in this approach. Consider, e.g., the intermediate state 6p = 5p 6 6p. It 
contains an admixture of states 5p 5 nsQd: Qp = 5p e 6p + a5p 5 ns6d + .... This mixed state is 
included into the sum (p7|). However, the sum (j2T|) must include all many-body states of 
opposite parity. This means that the state 5p 5 ns6d = 5p 5 ns6d — a5p e 6p + ... should also 
be included into the sum. Such contributions to E PNC have never been estimated directly 
within the sum-over- states approach. However, they are included into our mixed-states 
calculation. The result of the sum-over-states approach, 0.909, is very close to the result 
of the mixed-states approach, 0.908. It is important to note that the omitted higher-order 
many-body corrections are different in these two methods. This may be considered as an 
argument that the omitted many-body corrections in both calculations are small. Of course, 
here we assume that the omitted many-body corrections to both values (which, in principle, 
are completely different) do not "conspire" to give exactly the same magnitude. 

Therefore we will take 0.908 for the value of Ep^c (Subtotal of Table |) as this corre- 
sponds to the most complete mixed-states calculation and is in agreement with the sum- 
over-states calculation of reference ||. 
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With Breit, our result becomes 0.902 x 10 ll ieap.(— Qw/N). This correction is in agree- 
ment with [0-|l4 . 

We use the two-parameter Fermi model for the proton and neutron distributions: 

p(r) =p [l + exp[(r-c)/a]]~ 1 , (28) 

where t = a(41n3) is the skin-thickness, c is the half-density radius, and po is found from 
the normalization condition / pdV = 1. In 1989 the thickness and half-density radius for 
the proton distribution were taken to be t p = 2.5 fm and c p = 5.6149 fm (corresponding to 
a root-mean-square (rms) radius (r 2 ) 1 ^ 2 = 4.836 fm). In this work we have used improved 
parameters t p = 2.3 fm and c p = 5.6710 fm ({r 2 ) 1 ^ 2 = 4.804 fm) f2"6fl . This changes the 
wave functions slightly, leading to a very small correction to the PNC amplitude of 0.08% 
(0.0007). (This is in agreement with a simple analytical estimate: the factor accounting for 
the change in the electron density is ~ (4.804/4.836)-^ a ~ 0.1% .) This correction has 
already been included into the TDHF value. 

In the work f3| we used the proton distribution in the weak interaction Hamiltonian 
(Eq. fllTD). In the current work we have found the small correction to Ep^c which arises 
from taking the (poorly understood) neutron density in Eq. (|i~7D . We use the result of 
Ref. [E7J] for the difference Ar np = 0.13(4) fm in the root-mean-square radii of the neutrons 
( r n) 1 in d protons (r 2 ) 1 / 2 . We have considered three cases which correspond to the same 
value of (r 2 ): (i) c n = c p , a n > a p ; (ii) c n > c p , a n > a p ; and (iii) c n > c p , a n = a p (using 
the relation (r 2 ) « \c 2 n + f 7r 2 a 2 ). We have found that E PNC shifts from -0.18% to -0.21% 
when moving from the extreme c n = c p to the extreme a n = a p . Therefore, Ep^c changes 
by about —0.2% (—0.0018) due to consideration of the neutron distribution. This is in 



agreement with Derevianko's estimate, —0.19(8)% |19| . 

In the next section we discuss the radiative corrections to Epnc- The dominating con- 
tribution comes from the Uehling potential and increases the amplitude by 0.4%. 
Therefore, we have 



E PNC = 0.9041 x l0- u iea B (-Q w /N) (29) 
as our central point for the PNC amplitude. The error will be estimated in Section 0. 



IV. QED-TYPE RADIATIVE CORRECTIONS TO ENERGY LEVELS, WAVE 
FUNCTIONS, AND THE PNC AMPLITUDE 

The radiative corrections to the weak charge Qw have been calculated for the free elec- 
tron. However, an electron in a heavy atom is bound, and this produces additional radiative 
corrections proportional to a(Za) n , n = 1,2, .... Recently such corrections were considered 
by Milstein and Sushkov |17j . They found that the most important are corrections enhanced 
by the large parameter ln(A/i?), where A = h/mc is the electron Compton wavelength and 
R is the nuclear radius. This type of correction arises from the radiative corrections to 
the electron wave function near the nucleus. In this region the s-wave and pi/2-wave (lower 
Dirac component) electron densities are singular, |?/>(r)| 2 ~ r ~z 2 o? ^ r^^ e ra di a tive corrections 
modify the potential at small distances r < A, V(r) = —Za(l +5)/r. Correspondingly, the 
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electron wave functions change, |^(r)| 2 ~ r z a ( 1+<5 ) for r < A. This gives the radiative 
correction factor for the electron density inside the nucleus, 



\lp{R)\ 2 (\\Z 2 c?2& 



l^(A)| 2 



exp(25ZVln(A/i?)) . (30) 



For the Uehling (vacuum polarization) potential 5 ~ a ln(A/r) [28|. This gives an additional 
power of the large parameter ln(A/i?). This led Milstein and Sushkov [I]J to conclude that 
the Uehling potential gives the dominating radiative correction to E PNC , ~ Z 2 a 3 \n 2 (X/ R). 
Numerical calculations of the Uehling potential contribution have been performed in [TEj 
and in the present work. This radiative correction increases Ep^c by 0.4%. 

Milstein and Sushkov JlTj demonstrated that there are no other radiative corrections 
which are enhanced by In (X/R). However, any correction to the potential with nonzero 
S(R) ~ a gives a correction to the electron density ~ Z 2 a 3 \n(X/R). There are also correc- 
tions ~ Z 2 a 3 \n(Z 2 a 2 ) which originate from the shift of the energy levels (Lamb shift). We 
briefly discuss these corrections below. 

Let us start our discussion from the radiative corrections to energy levels. The calculation 
of the shift can be divided into two parts: one in which the electron interaction with virtual 
photons of high-frequency are considered, and one in which virtual photons of low-frequency 
are considered. 

In the high-frequency case the external field (the strong nuclear Coulomb field) need only 
be included to first-order. In this case the contributions to the Lamb shift arise from the 
diagrams presented in Fig. |[ The contribution of the Uehling potential (Fig. §(a)) to the 
Lamb shift is very small. The main contribution comes from the vertex correction (Fig. |8](b)). 
(In the case of a free electron the vertex diagrams give the electric f{q 2 ) and magnetic g(q 2 ) 
form factors.) The perturbation theory expression for f(q 2 ) contains an infra-red divergence 



and requires a low-frequency cut-off parameter k - see, e.g., [p8| . Assuming q 2 <C m 2 c 2 , the 



high-frequency contribution to the Lamb shift can be presented as a potential given by the 



following expression [28 
5$(r) 



5$/ + 5$u 
ah 2 r , m 11 1\ A T „ . . ah 



, J n— + )A$(r)-« 7-V$(r). (31) 

37rm 2 c 2 V 2k 24 5^ V; 4nmc 1 K) K J 

For the Coulomb potential, A$ = —4irZe5(r). Here the last long-range term (S$> g ) comes 
from the anomalous electron magnetic moment (g(0)); the infra-red cut-off parameter k 
appears from the electric form factor f(q 2 )- This term with the large In ^ gives the dominant 
contribution to the Lamb shift of s-levels. The infra-red divergence for k — > indicates the 
importance of the low-frequency contribution for this term. 

If we go beyond the approximation q 2 << m 2 c 2 , the term 5Qf should be associated with a 
non-local self-energy operator S ra d(r, r', E) with typical values |r — r'| < — and r ~ r' ^ 
We need this operator in a simple limit, E « mc 2 . Matrix elements of E ra d(r, r', 0) depend 
only on the electron density near the origin. Therefore, we can express the Lamb shift of the 
external electron state in a neutral atom in terms of the known Lamb shift of highly excited 
states in hydrogen- like ions. To implement this scheme we can approximate S ra d(r, r', 0) by 
a two-parametric (A, b) potential 
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5 $ = -A^re-v^rQ ; (32) 

7T 

where the factor r is introduced to remove the singularity at r = 0. (We have also performed 
the calculation using the potential 5<&f without this factor r. The results for the energy 
levels are the same.) The parameters A and b in 5Qf can be found from the fit of the Lamb- 
shift of the high Coulomb levels 3s, 4s, 5s and 3p, Ap and bp (in one-electron ions) which 
were calculated as a function of the nuclear charge Z in Refs. |[29|| . We have checked that 
A s = 180, A p = 90 and b = 1 fit all these Lamb shifts quite accurately (for the potential 
without the factor r we have found A s = 1.17 and A p = 1.33). The anomalous magnetic 
moment contribution is 

5$ = • V$ . (33) 

Airnic 

The Uehling potential for a finite nucleus is given by |30| (in atomic units: % — m = e = 1, 

a = 1/c) 

= -— / dx X P (X) / dt 4 + -±-)U- 2t \ r -*\/ a - e -2t(r+x)/a\ ,34) 

3r Jo 7i ^t d 2t 5/ v 7 

where p(x) is the nuclear charge density. It is more convenient to use a simpler formula for 
5&u for r > R, 

S<S>u(r) = $(r)^f(ft V^T(i + ^)e^/(x) , (35) 
/(x) = -e x + e~ x + xe x + xe~ x , x = 2tR/a , (36) 

and take ^y^p = ^ $(" r = • There is practically no loss of numerical accuracy in this 
approximation since a typical scale for the variation of 5&u/& is given by the electron 
Compton length — >> R. 

Since = 0.4 is not so small for Cs it is important to estimate the contribution of the 
higher-order in Za vacuum-polarization correction (the Wichmann-Kroll term [[H|]). For 
simplicity we use an approximate expression for this potential 

2 a 0.092Z 2 «2 

5<& WK = 1 flQ a ® (37) 

3vr 1 + (L62r) 4 v ; 

which is exact at small and large distances (a small-r asymptotic was presented in Refs. 



32|,[17| ). We have found that the contribution of this potential to the s-wave energies is 
~ 30 times smaller than that of the Uehling potential. The contribution of the Uehling 
potential 5§u to the Lamb shift is always small (less than 15%). Note that the contribution 
of the radiative corrections to the electron core electrostatic potential can be estimated 
using the semiclassical expression for E ra d(r, r', E) presented in Ref. [Q. This contribution 
for s-levels is two orders of magnitude smaller than that of the nuclear potential (for higher 
angular momenta the nuclear contribution is small and the electron contribution is relatively 
important). 

The potential <5$ 9 due to the magnetic form factor is a long-range one. This also hints 
that there are no large higher Za corrections here. The contribution of this potential to 
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the Lamb shift is about 30% for s-levels and 70% for p-levels. Note that for the term 5<&f 
and the total Lamb-shift we do not use the assumption Za « 1 since we fitted the exact 
results for the single-electron ions. 

The radiative corrections for Cs energy levels are presented in Table |I[ 

Now we discuss the contribution of QED-type radiative corrections to the electron wave 
function and PNC amplitude Ep^c- Note that it is not enough to calculate the radiative 
corrections to the matrix element of the weak interaction (n'pi/2\Hw\ ns ) ■ Corrections to 
the energy intervals like 6s — 6p are also important since these intervals are small at the 
scale of the atomic unit (~ 1/20) and sensitive to perturbations. The change in the energies 
also influences the large- distance behavior of the electron wave functions which determine 
the usual El amplitudes in the sum-over-states approach. If we keep only three dominating 
terms in the sum-over- states approach (see below) the contributions of the energy shifts 
(-0.29 %) and corrections to the amplitudes (0.33%) cancel each other. The simplest way to 
find the total answer (including the corrections to the weak matrix elements) is to include 
<5$ into the relativistic Hartree-Fock equations and then perform all calculations. 

The results are the following: the largest contribution to Ep^c comes from the Uehling 
potential 5&u- It increases Ep^c by 0.41% (in agreement with |T8|jr7[| ). The contribution of 
5§ g is —0.03% (due to cancellation of the contributions of the corrections to the s-wave and 
p-wave); the Wichmann-Kroll contribution is —0.006%. Note that the latter contribution 
can be estimated analytically using Eq. fl3~0l) . The ratio of the Wichmann-Kroll contribution 
to the Uehling contribution in the logarithmic approximation is —0.184 Z 2 a 2 /ln(X/R). 

The contribution of <5$/ is well-defined only when it originates from large distances (due 
to the shift of the energy levels and the large-distance behavior of the electronic wave function 
which influences electromagnetic amplitudes). A small-distance contribution is not gauge- 
invariant and should be treated simultaneously with the correction to the Z-boson exchange 



vertex [ 17J . Therefore, we deliberately selected a non-singular parametric potential ( p2|) 
to approximate 5$/. This potential does not generate any singular (in nuclear radius R) 
contributions proportional to the large parameter ln(A/i?). This means that the subject of 
our discussion now is different from the radiative corrections to the weak matrix element 
considered by Milstein and Sushkov in [|T7|]. They were interested in contributions enhanced 
by this large parameter hx(X/R). 

The main contribution to the Lamb shift comes from the distances r < 1/Z (in atomic 
units). This corresponds to the parameter b > Za in the potential (0). On the other 
hand a typical distance for the radiative corrections is r ~ A = — . Therefore, to estimate 
the contribution of 5$ f we performed calculations for two extreme values of the parameters 
6 = 1 and b = Za (the values of the Lamb shifts are kept the same in both cases by selecting 
appropriate values of the parameter A). For 6 = 1 this contribution is —0.2%. For b = Za 
it is -0.08%. This gives us an indication that the large-distance contribution to the radiative 



corrections to Ep^c is small. On the other hand, according to Sushkov and Milstein ||T7 
the small-distance contribution to E PNC is dominated by the Uehling potential which was 
considered above. 
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V. ESTIMATE OF ACCURACY OF PNC AMPLITUDE 



We have estimated the error of the PNC amplitude in a number of different ways. There 
are two main methods: (i) root-mean-square (rms) deviation of the calculated energy in- 
tervals, El amplitudes, and hyperfme structure constants from the accurate experimental 
values; (ii) influence of fitting of energies and hyperfme structure constants on the PNC 
amplitude. 



A. Root-mean-square deviation 

Remember that the PNC amplitude can be expressed as a sum over intermediate states 
(see formula (27)). Each term in the sum is a product of El transition amplitudes, weak 



matrix elements, and energy denominators. There are three dominating contributions to 
the 6s — 7s PNC amplitude in Cs: 

(7s\H E1 \6p)(6p\H w \6s) (7s\H w \6p) (6p\H E1 \6s) (7s\H E1 \7p) (7p\H w \6s) 

hpNC ~ W W + P W F F '•• 

= -1.908 + 1.493 + 1.352 + ... = 0.937 + ... (38) 

(the numbers are from the work || where the sum-over-states method was used; here we 
just demonstrate that these terms dominate). While we do not use the sum-over-states 
approach in our calculation of the PNC amplitude, it is instructive to analyze the accuracy 
of the El transition amplitudes, weak matrix elements, and energy intervals which contribute 
to Eq. ( [38]) as they have been calculated using the same method as that used to calculate 

EpNC- 

Let us begin with the energy intervals. The calculated removal energies are presented in 



Table |TJ. The Hartree-Fock values deviate from experiment by 10%. Including the second- 
order correlation corrections E^ reduces the error to ~ 1%. When screening and the hole- 
particle interaction are included into S*- 2 -* in all orders, the energies improve, ~ 0.3%. The 
percentage deviations from experiment of the energy intervals of interest are: Eq s — Eq p , 0.3; 
E 7s — Eq p , 0.4; and E 6s — E 7p , 0.3. The rms error is 0.3%. We can in fact reproduce energy 
intervals exactly by placing coefficients before the correlation potential. We will use this 
procedure as another test of the stability of the results. Note, however, that the accuracy 
for the energies is already very high and the remaining discrepancy with experiment is of 
the same order of magnitude as the Breit and radiative corrections. Therefore, generally 
speaking, we should not expect that fitting of the energy will always improve the results for 
amplitudes and hyperfme structure. In fact, as we will see below, some values do improve 
while others do not. The overall accuracy, however, remains at the same level. 

The relevant El transition amplitudes (radial integrals) are presented in Table [IV]. These 
were calculated with the energy-fitted "bare" correlation potential E^ 2 ) and the (unfitted and 
fitted) "dressed" potential E. Structural radiation and normalization contributions were also 
included. In Table [V] the percentage deviations of the calculated values from experiment are 
listed. Without energy fitting, the rms error is 0.1%. Fitting the energy gives a rms error 
of 0.2% for £W an d q.3% for the complete E. 
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We cannot directly compare weak matrix elements with experiment. Like the weak 
matrix elements, hyperfine structure is determined by the electron wave functions in the 
vicinity of the nucleus, and this is known very accurately. The hyperfine structure constants 
calculated in different approximations are presented in Table |VI[ Corrections due to the 
Breit interaction, structural radiation, and normalization are included. The percentage 
deviations from experiment are shown in Table [VI \ The rms deviation of the calculated 
hfs values from experiment using the unfitted E is 0.5%. With fitting, the rms error in the 
pure second-order approximation is 0.3%; with higher orders we get 0.4%. We are, however, 
trying to estimate the accuracy of the s — p weak matrix elements. It seems reasonable for 



us to use the square-root formula, y hfs(s)hfs(p). The errors are presented in Table |V11 . 
Notice that by using this approach the error is smaller. Without energy fitting, the rms 
error is 0.5%. With fitting, the rms error in the second-order calculation (S^) is 0.2% and 
in the full calculation (E) it is 0.3%. 

From this section we can conclude that the rms error for the relevant parameters is 0.5% 
or better. 

Note that from this analysis the error for the sum-over-states calculation of -EWc* would 
be larger than this, as the errors for the energies, hfs constants, and El amplitudes contribute 
to each of the three terms in Eq. (|3^). However, in the mixed-states approach, the errors 
do not add in this way. We get a better indication of the error of our calculation of Ep^c 
in the next section. 



B. Influence of fitting on the PNC amplitude 

In the section above we presented calculations in three different approximations: with 
unfitted E, and with E^ and E fitted with coefficients to reproduce experimental removal 
energies. The errors in these approximations are of different magnitudes and signs. We now 
calculate the PNC amplitude using these three approximations. The spread of the results 
can be used to estimate the error. 



The results are listed in Table |X[ It can be seen that the PNC amplitude is very 
stable. The PNC amplitude is much more stable than hyperfine structure. This can be 
explained by the much smaller correlation corrections to Ep^c (~ 2% for Ep^c and ~ 30% 



for hfs; compare Table | with Table |VlD . One can say that this small value of the correlation 



correction is a result of cancellation of different terms in (^) but each term is not small 
(see Table |). However, this cancellation has a regular behavior. The same correlation 
potential E is used to calculate energies and correlation corrections ( |2T)| ) to E PNC . Therefore, 
whatever way is used to calculate E, if the accuracy for energies is good, the correlation 
correction is stable. The stability of Ep^c m &y be compared to the stability of the 
usual electromagnetic amplitudes where the error is very small (even without fitting). 

We have also considered the fitting of hyperfine structure using different coefficients 
before each E. The first-order in E correlation correction fl26| ) changes by about 10%. This 
changes the PNC amplitude by about 0.4%. 

It is also instructive to look at the spread of Ep^c obtained in different schemes. (This 



has already been discussed in some detail in Section |TJ.) The result of the present work is 



in excellent agreement with our earlier result 0] while the calculation scheme is significantly 
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different. The only other calculation of the Ep^c hi Cs which is as complete as ours is that of 
Blundell et al. [||]. Their result in the all-orders sum-over-states approach is 0.909 (without 
Breit) and is very close to our value of 0.908 (corresponding to "Subtotal" of Table |). Our 
result (0.904) obtained in second-order with fitting of the energies is useful in determining 
the accuracy of the calculations of Ep^c- (Remember that this value is in agreement with 
results of similar calculations performed in PJl^l; see Section |T|.) One can see that replacing 
the all-order £ by its very rough second-order (with fitting) approximation changes Ep^c 
by less than 0.4% only. On the other hand, if the higher orders are included accurately, the 
difference between the two very different approaches is 0.1% only. 

The maximum deviation we have obtained in the above analysis is 0.5%. We will therefore 
use this as our estimate for the uncertainty of the Ep^c calculation. 

We do not include the error associated with the radiative corrections into the estimate 
of the accuracy since the large-distance (r > A) contribution of the radiative corrections 
is small and the small-distance contribution (R < r < A) enhanced by ln(A/i?) can be 
accurately calculated. 



VI. CONCLUSION 



We have obtained the result 



E 



PNC 



0.904 1 ± 0.5% x l0- ll iea B (-Q w /N) 



(39) 



for our calculation of the 6s — 7s PNC amplitude in Cs. This is in agreement with other PNC 
calculations, however we would like to emphasize that our calculation is the most complete. 
The most precise measurement of the 6s — 7s PNC amplitude in Cs is B 



lm(E 



PNC 







1.5939(56) 



mV 



cm 



(40) 



where (3 is the vector transition polarizability. There are currently two very precise values 
for p. One value 



was obtained in our analysis 
M1 M J(3 ratio 1D| 



P = 26.957(51)4 

of the Bennett and Wieman measurements 
We have obtained another value 



P = 27.15(ll)a| 



(41) 



of the 



(42) 



from the measurement ^1] oia/P and an analysis, similar to that of Ref. |]2*I |, using the most 
accurate experimental data for the El transition amplitudes including recent measurements 
of Vasilyev et al. ||42| . The errors in Eqs. (1"I,4"2"D are obtained by adding in quadrature the 



experimental and theoretical errors. Notice that the central point of our value (Eq. (|42]) ) 
differs slightly from the value 27.22(11) obtained in the work | 4"2|| . Note also that the value 



(f!2|) coincides with the value presented in |2T[ which was obtained with slightly different El 
amplitudes. This is because of an accidental cancellation of the changes to different terms. 
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Using the conversion \e\/a 2 B = 5.1422 x 10 12 mV/cm, we therefore obtain for the weak 
charge of the Cs nucleus with j3 = 26.957: 

Q w = -72.09(29) cxp (36) theor , (43) 

or with f3 = 27.15 

Q w = -72.60(39) cxp (36) theor , (44) 

where the experimental error is obtained by adding in quadrature the error for (5 and the 
error for Im(E PNC )/ ft. 

If we take an average value for /3 

(3 = 26.99(5)4 ( 45 ) 

then 



Q 



w 



-72.18(29) exp (36) theor 



(46) 



These results (Eqs. ( ^3| , ^4| , ^6|) ) deviate by 2.2a, 0.9a and 2.0a, respectively, from the 
standard model value Qw = —73.09(3) |22| (note that the standard deviations a are different 
in each case). 



ACKNOWLEDGMENTS 



We are grateful to A. Milstein, O. Sushkov, and M. Kuchiev for useful discussions. This 
work was supported by the Australian Research Council. 



17 



REFERENCES 



[1] LB. Khriplovich, Parity Nonconservation in Atomic Phenomena (Gordon and Breach, 
Philadelphia, 1991). 

M.-A. Bouchiat and C. Bouchiat, Rep. Prog. Phys. 60, 1351 (1997). 

C. S. Wood et al. } Science 275, 1759 (1997). 

V.A. Dzuba, V.V. Flambaum, and O.P Sushkov, Phys. Lett. A 141, 147 (1989). 
S.A. Blundell, W.R. Johnson, and J. Sapirstein, Phys. Rev. Lett. 65, 1411 (1990); S.A. 
Blundell, J. Sapirstein, and W.R. Johnson, Phys. Rev. D 45, 1602 (1992). 
R.J. Rafac, and C.E. Tanner, Phys. Rev. A, 58 1087 (1998); R.J. Rafac, C.E. Tanner, 
A.E. Livingston, and H.G. Berry, Phys. Rev. A, 60 3648 (1999). 
S.C. Bennett, J.L. Roberts, and C.E. Wieman, Phys. Rev. A 59, R16 (1999). 
S.C. Bennett and C.E. Wieman, Phys. Rev. Lett. 82, 2484 (1999); 82, 4153(E) (1999); 
83, 889(E) (1999). 

R. Casalbuoni, S. De Curtis, D. Dominici, and R. Gatto, Phys. Lett. B 460, 135 (1999). 
J. L. Rosner, Phys. Rev. D 61, 016006 (1999). 
J. Erler and P. Langacker, Phys. Rev. Letts. 84, 212 (2000). 
A. Derevianko, Phys. Rev. Lett. 85, 1618 (2000). 

V.A. Dzuba, C. Harabati, W.R. Johnson, and M.S. Safronova, Phys. Rev. A 63, 044103 
(2001). 

M.G. Kozlov, S.C Porsev, and I.I. Tupitsyn, Phys. Rev. Lett. 86, 3260 (2001). 
W.J. Marciano and A. Sirlin, Phys. Rev. D 27, 552 (1983); W.J. Marciano and J.L. 
Rosner, Phys. Rev. Lett. 65, 2963 (1990). 

BW. Lynn and P.G.H. Sandars, J. Phys. B. 27, 1469 (1994); I. Bednyakov et al, Phys. 
Rev. A 61, 012103 (1999). 

A.I. Milstein and O.P. Sushkov, e-print |hep-ph/0109257| . 

W.R. Johnson, I. Bednyakov, and G. SofT, Phys. Rev. Lett. 87, 233001 (2001). 
A. Derevianko, Phys. Rev. A 65, 012106 (2002). 
V.A. Dzuba and V.V. Flambaum, Phys. Rev. A 62, 052101 (2000). 
V.A. Dzuba, V.V. Flambaum, and O.P. Sushkov, Phys. Rev. A 56, R4357 (1997). 

D. E. Groom et al, Euro. Phys. J. C 15, 1 (2000). 

V.A. Dzuba, V.V. Flambaum, PC Silvestrov, and O.P. Sushkov, J. Phys. B 20, 1399 
(1987). 

V.A. Dzuba, V.V. Flambaum, and O.P. Sushkov, Phys. Lett. A 140, 493 (1989). 
V.A. Dzuba, V.V. Flambaum, A.Ya. Kraftmakher, and O.P. Sushkov, Phys. Lett. A 
142, 373 (1989). 

G. Fricke et al, At. Data and Nucl. Data Tables 60, 177 (1995). 
A. Trzcihska et al, Phys. Rev. Lett. 87, 082501 (2001). 

V.B. Berestetskii, E.M. Lifshitz, and L.P. Pitaevskii, Relativistic Quantum Theory 
(Pergamon Press, Oxford, 1982). 

P.J. Mohr and Y.-K. Kim, Phys. Rev. A 45, 2727 (1992); P.J. Mohr, Phys. Rev. A 46, 
4421 (1992). 

LW. Fullerton and G.A. Rinker, Jr., Phys. Rev. A 13, 1283 (1976). 

E. H. Wichmann and N.M. Kroll, Phys. Rev. 101, 343 (1956). 
A.I. Milstein and V.M. Strakhovenko, ZhETF 84, 1247 (1983). 
V.V. Flambaum and V.C Zelevinsky, Phys. Rev. Lett. 83, 3108 (1999). 



18 



[34] C.E. Moore, Natl. Stand. Ref. Data Ser. (U.S., Natl. Bur. Stand.), 3 (1971). 
[35] R.J. Rafac, C.E. Tanner, A.E. Livingston, and H.G. Berry, Phys. Rev. A 60, 3648 
(1999). 

[36] M.-A. Bouchiat, J. Guena, and L. Pottier, J. Phys. (France) Lett. 45, L523 (1984). 

[37] E. Arimondo, M. Inguscio, and P. Violino, Rev. Mod. Phys. 49, 31 (1977). 

[38] S.L. Gilbert, R.N. Watts, and C.E. Wieman, Phys. Rev. A 27, 581 (1983). 

[39] R.J. Rafac and C.E. Tanner, Phys. Rev. A 56, 1027 (1997). 

[40] M.-A. Bouchiat and J. Guena, J. Phys. (France) 49, 2037 (1988). 

[41] D. Cho et al, Phys. Rev. A 55, 1007 (1997). 

[42] A. A. Vasilyev, I.M. Savukov, M.S. Safronova, and H.G. Berry, e-print jphysics/0112071 . 



19 



TABLES 



TABLE I. Contributions to the 6s — 7s Ep^c amplitude for Cs in units 10 ll ieaB{—Qw/N). 
(S corresponds to the (unfitted) "dressed" self-energy operator.) 



TDHF 0.8898 

(?Jj 7s \± s (e 7s )\6X 6s ) 0.0773 

(5ip 7s \t p (e 7s )\X 6s ) 0.1799 

{5Y 7s \t s (e 6s )^ 6s ) -0.0810 

(Y 7s \t p (e 6s )\5f 6s ) -0.1369 

Nonlinear in £ correction -0.0214 

Weak correlation potential 0.0038 

Structural radiation 0.0029 

Normalization -0.0066 

Subtotal 0.9078 

Breit -0.0055 

Neutron distribution -0.0018 

Radiative corrections 0.0036 

Total 0.9041 



TABLE II. Radiative corrections to RHF removal energies; units cm L (See also Table [Hi] .) 

_6s 7s 6pi/ 2 7pi/ 2 

-18.4 -5.0 0.88 0.31 



TABLE III. Removal energies for Cs in units cm . 



State 



RHF 



S(2) 



Experiment 



6s 27954 

7s 12112 

6p 1/2 18790 

7pi/ 2 9223 

a Taken from [34]. 



32415 
13070 
20539 
9731 



31492 
12893 
20280 
9663 



31407 
12871 
20228 
9641 
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TABLE IV. Radial integrals of El transition amplitudes for Cs in different approximations. 
The experimental values are listed in the last column, (a.u.) 



Transition 


RHF 


TDHF £W 

with fitting 


S 


t 

with fitting 


Experiment 


6s — 6p 
7 s — Qp 
7s — 7p 


6.464 
5.405 
13.483 


6.093 5.499 
5.450 5.198 
13.376 12.602 


5.497 
5.190 
12.601 


5.509 
5.204 
12.612 


5.497(8) a 
5.185(27) b 
12.625(18) c 


a Ref. g|. 
b Ref. [36]. 
c Ref. @. 












TABLE V. Percental 
different approximations. 


^e deviation from experiment of calculated El transition amplitudes in 


Transition 




S( 2 ) 

with fitting 


Percental 


$e deviation 
t 


t 

with fitting 


6s — 6p 
7 s — 6p 
7s — 7p 




0.04 
0.3 
-0.2 




0.0 
0.1 
-0.2 


0.2 
0.4 
-0.1 



TABLE VI. Calculations of the hyperfine structure of Cs in different approximations. In the 
last column the experimental values are listed. Units: MHz. 



State 


RHF 


TDHF 


S(2) 

with fitting 


t 


S 

with fitting 


Experiment 


6s 


1425.0 


1717.5 


2306.9 


2315.0 


2300.3 


2298.2 a 


7s 


391.6 


471.1 


544.4 


545.3 


543.8 


545.90(9) b 


6P1/2 


160.9 


200.3 


291.5 


293.6 


290.5 


291.89(8) c 


7P1/2 


57.6 


71.2 


94.3 


94.8 


94.1 


94.35 a 



a Ref. @. 
b Ref. |8|. 
c Ref. |§. 
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TABLE VII. Percentage deviation from experiment of calculated hyperfine structure constants 
in different approximations. 



State 


S( 2 ) 

with fitting 


Percentage deviation 
E 


E 

with fitting 


6s 


0.4 


0.7 


0.09 


7s 


-0.3 


-0.1 


-0.4 


6pi/ 2 


-0.1 


0.6 


-0.5 


7pi/ 2 


-0.05 


0.5 


-0.3 



TABLE VIII. Percentage deviation from experiment of calculated v /hfs(s)hfs(p) (we will denote 



this by s — p in the tables) in different approximations. 



v /hfs(s)hfs(p) 


S( 2 ) 

with fitting 


Percentage deviation 
E 


S 

with fitting 


6s — 6p 


0.1 


0.7 


-0.2 


6s - 7p 


0.2 


0.6 


-0.09 


7s — 6p 


-0.2 


0.2 


-0.4 



TABLE IX. Values for Ep^c m different approximations; units 10 n ieaB(— Qw/N). 





£( 2 ) with fitting 


S 


E with fitting 


Epnc 


0.901 


0.904 


0.903 
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(a) (b) 
FIG. 1. Second-order (a) direct and (b) exchange correlation diagrams in the Feynman diagram 
technique. Dashed line is the Coulomb interaction. Loop is the polarization of the atomic core. 



+ .... 



FIG. 2. Hole-particle interaction in the polarization operator. 



FIG. 3. Screening of the Coulomb interaction in the direct correlation diagram. 



FIG. 4. The screened Coulomb interaction (with hole-particle interaction). 
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FIG. 5. The electron self-energy operator with screening and hole-particle interaction included. 
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FIG. 6. Lowest-order correlation corrections to the PNC El transition amplitude. Dashed line 
is the El field; cross is the nuclear weak field. 
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FIG. 7. Small corrections to the PNC El transition amplitude: external field inside the correla- 
tion potential. In diagrams (a) the weak interaction is inside the correlation potential (ST, denotes 
the change in £ due to the weak interaction); this is known as the weak correlation potential. 
Diagrams (b) represent structural radiation (photon field inside the correlation potential). 




(a) (b) 

FIG. 8. High-frequency contribution to radiative corrections. Diagram (a) corresponds to the 
Uehling potential. Diagram (b) is the vertex correction. The single solid line represents the bound 
electron; the double line is the free electron; the Coulomb interaction is denoted by the dashed 
line; and the filled circle denotes the nucleus. 
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